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ABSTRACT 

We report the discovery of 40 bright quasars and active galactic nuclei (AGNs) 
at low Galactic latitude (b < 1 20° | ) . The low galactic latitude region has been 
considered as a place to avoid when searching for extragalactic sources, because 
of the high Galactic extinction as well as a large number of stars contaminating 
the sample selection. Bright quasars (R < 17) suffer more from such difficulties 
because they look like bright stars, which are numerous at low b, yet their surface 
number density is very low. In order to find quasars in this region of the sky 
less explored for extragalactic sources, we have started a survey of low Galactic 
latitude bright quasars as a part of the Seoul National University Quasar Sur- 
vey in Optical (SNUQSO). Quasar candidates have been selected from radio and 
near-infrared (NIR) data. Out of 88 targets, we identify 29 bright quasars/AGNs 
around the antigalactic center, and 11 bright quasars/AGNs in the outskirts of 
the Galactic center, from two observing runs in 2006 at the Bohyunsan Opti- 
cal Astronomical Observatory (BOAO) in Korea Our finding demonstrates that 
quasars/AGNs can be discovered effectively even at low Galactic latitude using 
multi-wavelength data. 

Subject headings: quasars: general — galaxies: active — techniques: spectro- 
scopic — surveys 
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1. INTRODUCTION 

Since the discovery of quasar 3C273 in 1963 (Schmidt 1963), more than 100,000 quasars 
and active galactic nuclei have been found to date (Schneider et al. 2005; Croom et al. 
2004; Hagen et al. 1995; Hewett et al. 1995; Becker et al. 2001; Schmidt & Green 1983; 
Veron-Cetty & Veron 2006). Despite of the rapid increase in the number of known quasars, 
however, only a small number of quasars have been discovered near the Galactic plane. The 
dots and crosses in Figure. 1 represent the positions of quasars discovered to date, taken 
from Veron-Cetty & Veron (2006). Clearly, there is a lack of quasars/AGNs at low Galactic 
latitude (b < |20|). There are two reasons for this paucity of quasars at low Galactic latitude. 
First, there are simply too many stars that look like quasars at low Galactic latitude so that 
it is much more challenging to select quasar candidates in such a region than at a high 
Galactic latitude region. High Galactic extinction in the Galactic plane also causes a trouble 
when searching for quasars. Not only does Galactic extinction dim the light from quasars 
by several magnitudes or more, but it also affects the optical colors of quasars, making it 
difficult to isolate quasar candidates with otherwise highly efficient quasar selection methods 
employing optical color-color diagrams (e.g., Richards et al. 2004). 

Among quasars, bright quasars (V < 17 mag) get special attention. Because of their 
brightness, they are easy to use to investigate the intergalactic medium (e.g., Weymann et al. 
1981; Rauch 1998; Fan et al. 2006), the properties of supermassive blackholes and their host 
galaxies (e.g., Kaspi et al. 2000; Kim et al. 2006), and the physical state and distribution of 
gas in our Galaxy (Sternberg 2003; Sembach et al. 2003; Savage & Sembach 1996). Because 
of their usefulness, bright quasars have been searched for extensively in high Galactic latitude 
regions, especially in the northern hemisphere (Schmidt & Green 1983; Hagen et al. 1995; 
Hewett et al. 1995; Becker et al. 2001; Schneider et al. 2005). These surveys appear to 
have found the majority of bright quasars in their respective search areas, although a small 
number of bright quasars are still waiting to be discovered (Lee et al. 2008). 

Bright quasars at low Galactic latitude also have many astrophysical applications. They 
can allow us to study the gaseous component of matter in our galaxy from a unique point 
of view through their absorption lines, and offer an opportunity to study their host galaxies 
using bright stars close to them as guide stars for adaptive optics. They can also serve 
as reference points to examine the motion of stars in our Galaxy. However, finding bright 
quasars at low Galactic latitude seems to be a daunting job considering the aforementioned 
difficulties. A simple quantitative argument is given below to demonstrate how difficult it 
is to identify bright quasars at low Galactic latitude. The surface number density of bright 
quasars at V < 17 mag is 0.1 deg~ 2 , while the number density of stars at low Galactic 
latitude ranges from a few thousand (anti-galactic center) to a few tens of thousand per 
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square degree at around I = 0° and b = 10° (Robin et al.2003). This is 10 to 100 times 
more than the number density of stars at high Galactic latitude (Spagna et al. 1996). These 
numbers become even worse if we take into account Galactic extinction which dims the 
light from extragalactic objects by more than a magnitude in most areas at low Galactic 
latitude. Assuming 1 magnitude extinction, the apparent magnitude cut of V < 17 mag 
at low Galactic latitude would actually sample quasars with V < 16 mag if they were at 
high Galactic latitude. At V < 17 mag, the number density of quasars is only 0.01 deg~ 2 . 
Therefore, the quasar-to-star ratio of an apparent magnitude-limited sample at V < 17 
mag becomes 1 in a few hundred thousand to a few million, or 100 to 1000 times worse 
than at high Galactic latitude. At V < 17, the optical multi-color selection gives a quasar 
identification efficiency of about 90%, but at low Galactic latitude this efficiency would drop 
to 9% to 0.9% following the above argument. Clearly an alternative method is needed for 
selecting bright quasars at low Galactic latitude. 

In order to broaden our knowledge of bright quasars toward the zone of avoidance, we 
have started a survey of bright quasars at low Galactic latitude regions, as a part of the Seoul 
National University Quasar Survey in Optical (SNUQSO). Here, we report the discovery of 
40 bright quasars/AGNs at low Galactic latitude, using a selection algorithm designed to 
efficiently sample low Galactic latitude quasars. Throughout the paper, the magnitudes are 
Vega-based, unless otherwise mentioned. Also, we adopt the cosmological parameters of 
Vl m = 0.3, A = 0.7, and H = 70 km s" 1 Mpc' 1 . 



2. SELECTION OF CANDIDATES 
2.1. Method 

The key to successfully finding low Galactic latitude quasars lies in an efficient algorithm 
that can distinguish quasars from the plethora of stars and dust in the Galactic plane. To 
achieve such a goal, we use a combination of radio and near-infrared (NIR) data. For the 
radio data, we use the NRAO VLA Sky Survey (NVSS; Condon et al. 1998) data set, and 
for the NIR data, we use the Two Micron All Sky Survey (2MASS) data set (Skrutskie et 
al. 2006). These data sets are chosen because they provide nearly complete coverage of the 
Galactic plane to the depth sufficient for finding a fair number of bright quasars. The NVSS 
data covers the whole sky above 5 = —40° to the depth of 2.5 mJy in 21cm. The 2MASS 
data covers the entire sky to the depth of K = 15 mag. The reason for the use of the radio 
and NIR data sets is elaborated below. 

Our selection of candidates starts with identifying radio sources at \b\ < 20. As the name 
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"quasar" originates from the word "quasi-stellar radio star," it has been widely recognized 
that radio emitting point sources are very likely to be quasars or AGNs. The usual procedure 
for finding quasars using radio data is to find an optical counterpart to the radio source and 
perform spectroscopic observation to determine whether it is a quasar. In high Galactic 
latitude regions, it is not very difficult to identify the optical counterpart of a radio source, 
given that the radio data have a reasonably good positional accuracy. For example, the 
FIRST survey data show that the matching radius of 1" is sufficiently good for finding an 
optical counterpart (White et al. 2000). However, that approach is not directly applicable in 
this case. The beam size of the NVSS data is large (FWHM = 45") compared to FIRST, and 
the quoted positional accuracy of point sources varies from rms ~ 1" (for sources brighter 
than 15 mJy) to 7 " (at the survey limit of 2.5 mJy), forcing us to adopt a conservative 
2-er radius of 15" as a search radius to identify the optical counterpart. However, within a 
15" radius circle of a NVSS radio source at low Galactic latitude, there are on average ~6 
point sources down to the optical magnitude of R < 19 mag and sometimes more. 

In order to improve the quasar selection efficiency by weeding out stars near the radio 
position, we add a NIR selection method. Low redshift quasars are known to have red J — K 
colors due to a steeply rising NIR spectra (e.g., Glikman et al. 2006), and NIR colors have 
been used previously to select quasars and AGNs (e.g., Figure 5 in Barkhouse & Hall 2001; 
see also, Nelson et al. 1998; Cutri et al. 2000; Maddox & Hewett 2006). Therefore, we 
applied a NIR color cut on point sources near the radio position and selected objects with 
J — K > 1.4 as quasar candidates. The NIR cut of J — K > 1.4 has been found to be 
the optimal cut for minimizing the contamination from stars and galaxies by performing the 
above selection method on known quasars/galaxies/stars in the fourth data release of the 
Sloan Digital Sky Survey (SDSS; Adelman-McCarthy, et al. 2006). In general, the selected 
candidates have K < 14 mag and an optical magnitude of R < 19. With the NIR cut, the 
contamination rate of stars/galaxies on the SDSS sample with the above method is found 
to be on the order of a few percent or less, but at the expense of preferentially selecting low 
redshift quasars at z < 0.4 (see also, Chiu et al. 2007). A more complete description of the 
sample selection and its efficiency will appear in I. Lee et al. (2008, in preparation). 

Although our method can be applied to a very crowded region in the Galactic plane, we 
have limited the initial survey area to the following region so that we get a reasonably high 
quasar identification efficiency of ~ 50%. 



\b\ < 20° 
5° < \b\ < 20° 
10° < 161 < 20' 



at 120° < I < 240° 
at 40° < I < 120° & 240° < / < 320' 
at 0° < I < 40° & 320° < I < 360° 
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Using this selection method, we identify 601 candidates in the above area covering 
10919.9 deg 2 , 88 of which were observed over ~ 2700 deg 2 during 2006, as described in the 
next section. 



3. OBSERVATION 

For spectroscopic confirmation of quasars, we observed the quasar candidates using 
the longslit low-resolution spectrograph mounted on the Bohyunsan Observatory Eschelle 
Spectrograph (BOES) of the 1.8m telescope at Bohyunsan Optical Astronomical Observatory 
(BOAO; Kim et al. 2003). The spectroscopic observations were carried out on 2006 June 
17-24 and 2006 December 20 - 25, among which about four nights were in observable 
condition under the seeing of 2 - 3". Using a 150 per mm grating and a slit width of 3.6", 
we covered the wavelength range of 4200 - 8400 A at the spectral resolution of X/5 A ~ 
366, or 5v ~ 820 kms^ 1 . The observational setting is basically identical to the default 
observational setting of SNUQSO, which is described in more detail in Lee et al. (2008). 
The June run targeted bright quasar candidates in the outskirts of the Galactic center, while 
in the December run, we observed targets near the anti-galactic center. We observed 33 
candidates in June, and 55 candidates in December. The observed data were reduced with 
packages, following a standard procedure including flux calibration. We estimate the 
accuracy of the relative flux calibration (as a function of wavelength) to be roughly 15% (Lee 
et al. 2007a). We originally took optical magnitudes from the USNO-A2.0 catalog (Monet 
et al. 1998) for the quasar candidates, but some of the apparent magnitudes appeared to be 
too bright. In order to check the USNO magnitudes, we performed the multi-band imaging 
observation of objects that are identified as quasars/AGNs using the 1.5m telescope at the 
Mt. Maidanak Observatory in Uzbekistan during 2006 August 1-12 and also in 2006 
June/ July (M. Im et al. 2008, in preparation). The median seeing during the imaging run 
hovered around 0.7"to 1.0", and the observations were carried out in photometric condition. 
Our imaging observation reveals that the USNO magnitudes are too bright by more than 1 
mag for the four objects with USNO magnitude R < 15 (object numbers, 33, 36, 37, and 39 in 
Table 1), and we have revised the i?-band and U-band magnitude accordingly. In other cases, 
we find that the average difference in .R-band magnitude is i?(USNO) - i?(Maidanak) ~ 0.12 
mag. Table 1 lists the apparent magnitudes of the low Galactic latitude quasars/AGNs. The 
i?-band and 5-band magnitude of objects 29-39 in Table 1 are from our imaging observation, 



1 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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while the rest come from the USNO catalog. For some of the problematic cases, the blending 
of objects might be responsible for the bright USNO magnitude since we find neighbors 
within ~6"of the quasar. For the rest, it is not clear why the USNO magnitudes are so 
bright. Monitoring of the quasars/AGNs may reveal the temporal variation of brightness of 
such quasars/AGNs. Based on our imaging data of 11 objects, we note that the magnitudes 
in Table 2 originating from the USNO catalog may be too bright when R < 15 mag. 



4. RESULTS AND DISCUSSIONS 

From the observed spectra, quasars were identified as objects with redshifted, broad 
Balmer emission lines with their FWHM line widths exceeding 1000 kms~ l (White et al. 
2000; Schneider et al. 2003). With such a criterion, we find 40 quasars among 88 candidates. 
Application of an absolute magnitude cut of < —22 AB magnitude (Schneider et al. 2003) 
reduces this number to 370 

Figure 1 shows the distribution of the 40 new quasars/AGNs in Galactic coordinates 
(triangles and squares). Also plotted in the figure are quasars and active galaxies listed in 
Veron-Cetty & Veron (2006). We also plot quasars brighter than V = 17 mag as crosses. 
The lack of bright quasars at low Galactic latitude is apparent from this figure. Especially 
in the Galactic plane near the anti-galactic center (120° < I < 240° and \b\ < 8°), only 7 
bright quasars/AGNs have been reported previously. We more than double this number now 
by discovering 13 additional quasars/AGNs in this region. 

In Figure 2, we present the spectra of the 40 quasars/AGNs, while Table 1 summarizes 
the basic observed properties of these quasars. The extinction corrected absolute i?-band 
magnitude, M R , is calculated with the formula, 

M R = m R -5x log(d L (z)) - 25 - K R {z) - A R . 

Here, m R is the i?-band apparent magnitude before extinction correction (col. [8] in Ta- 
ble 1), d,L is the luminosity distance, K R (z) is the i^-correctionj^l derived from the composite 
spectrum of SDSS quasars (Vanden Berk et al. 2001), and A R is the Galactic extinction 



2 Mi can be calculated from the i?-band absolute magnitude listed in Table 1, by applying the color 
correction of i?(Vega) - %ab — -0.1 mag, which we derived using a template quasar spectrum from Vanden 
Berk et al. (2001). Here, i?(Vega) - iab would be ~ -0.2 mag for the i?-band magnitudes from USNO based 
on the comparison of the USNO i?-band magnitude to our data. 

3 Note that the correction is calculated assuming the Bessel i?-band filter response on the 4k x 4k CCD 
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toward the position of the quasar/ AGN taken from NASA/IPAC Extragalactic Database, 
which is based on the reddening map of Schlegel et al. (1998). Figure 2 shows a clear de- 
tection of redshifted Ha, B.f3, and [OIII] lines. Redshifts of these quasars/AGNs range from 
z = 0.05 to 0.3, and the rest-frame FWHMs are all measured to be well above 1000 kms^ 1 , 
ranging from 1300 to 8000 /cms -1 , showing that they belong in the category of quasars and 
active galaxies (White et al. 2000; Schneider et al. 2003). We also find that 13 are X-ray 
sources and 27 are infrared sources. In particular, SNUQSO J210931. 9-353258 has been 
previously identified as a gamma-ray blazar source by Sowards-Emmerd et al. (2005), but 
without redshift identification. We now identify this object to be a quasar at z — 0.202. 

The E(B — V) values in the direction of quasars/AGNs range from 0.1 to 0.9 with the 
median value at E(B — V) = 0.214. We examined whether there is a correlation between 
the Balmer decrement and the E(B — V) value. To do so, spectra are arranged according 
to the order of increasing E(B — V) value in Figure 2. The figure shows that there is a 
general tendency of weakening of Hf3 flux with respect to Ha flux as the E(B — V) values 
increase. We have plotted the Hf3 to Ha ratio versus E(B — V) and confirmed this general 
trend, but find a large scatter in the relation too. Detailed model-fitting of spectra and the 
observation during photometric nights may reduce the scatter, but it may not help much, 
since the large scatter could indicate that the line ratios of some of the quasars/AGNs do 
not reflect the simple case B recombination (e.g., Glickman et al. 2006; Soifer et al. 2004). 
The intrinsic scatter in this ratio due to dust absorption within quasars/AGNs and their 
host galaxies also complicates the interpretation. 

The radio-to-optical ratios (R*) of the newly discovered quasars/AGNs range between 

0. 1 and 282 with the median value at 1.9. The radio-to-optical flux ratio is defined to be the 
ratio of flux at the rest-frame 6 cm to 0.44 /iin (Stocke et al. 1992). The 6cm flux is derived 
from the 21cm flux assuming a power-law of f v ~ z/-°- 46 (Ho & Ulvestad 2001), and applying 
a corresponding power-law if-correction in the form of K{z) = —2.5 (a + 1) log{l + z) where 
a = —0.46. The 0.44 fim flux is derived from the i?-band apparent magnitude in Table 

1, corrected for Galactic extinction (Schlegel et al. 1998), as well as for the .^-correct ion 
derived from the SDSS quasar composite spectrum of Vanden Berk et al. (2001). We find 
six radio-loud quasars (R* > 10) which corresponds to 16% of our sample, excluding three 
AGNs. This number is slightly higher than the often-quoted value of 10% from more complete 
surveys at the magnitude range of our sample (e.g., Figure 18 in Ivezic et al. 2002; Hewett 
et al. 2001). The difference appears to be at least partially related to the radio selection of 



Camera installed on the 1.5m telescope at the Maidanak Observatory (Ko et al. 2007, in preparation). 
Application of this if-correction to the photographic i?-band magnitude from USNO may result in a difference 
on the order of ^0.1 mag. 
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our sample. 

In Figure 3, we plot the i?-band magnitude distribution of the newly discovered quasars/AGNs, 
before and after applying the extinction correction assuming the reddening map of Schlegel 
et al. (1998). Before applying the extinction correction, the distribution peaks at R ~ 16 
mag, with a wide tail toward R ~ 18 mag, showing that our sample completeness is skewed 
toward the brighter magnitude of it! < 16 mag. After applying the extinction correction, 
many R > 16 mag quasars/AGNs move toward the R < 16 mag area, making our sample ap- 
pear as if it were selected to be R < 16 mag. Note that the Galactic extinctions, Ar and Ab, 
range from 0.22 to 2.57 mag and 0.25 to 4.04 mag, respectively, toward the quasars/AGNs. 

Some of the quasars are found to be fairly bright. If USNO apparent magnitudes are ac- 
curate (see Section 3), SNUQSO J035657.03+425540.6 and SNUQSO J062355.21+001843.4 
would have R < 14 after Galactic extinction correction, making them two of the brightest 
quasars discovered to date. 

The 40 quasars/AGNs are discovered over about 2700 deg 2 of the sky, and therefore 
their surface number density is 0.014 deg~ 2 down to R < 16 mag. We compare this number 
to the surface density of optically selected quasars/AGNs at high Galactic latitude from La 
Franca & Cristiani (1997), and also to the radio selected quasars from White et al. (2000). 
Assuming B — R ~ 0.7 mag for quasars (Becker et al. 2001), our effective magnitude cut 
of R < 16 mag can be converted to B < 17 mag. At that magnitude limit, La Franca & 
Cristiani (1997) find the surface number density of quasars to be roughly 0.05 deg~ 2 (See also 
Croom et al. 2004). On the other hand, White et al. (2000) find about 70 quasars at R < 16 
over 2682 deg 2 which gives the surface density of 0.026 deg~ 2 . Our surface number density is 
nearly one half of the radio-selected sample of White et al. (2000). The major reason for the 
discrepancy is our preferential selection of lower redshift quasars/AGNs with the adopted 
J — K cut, which leads to the exclusion of about 50% of R < 16 quasars at z > 0.4 (see 
Figure 8 in White et al. 2000). Also, a part of the discrepancy can be explained with the 
difference in the depth of the radio data (2.5 mJy versus our 1 mJy). The discrepancy in 
the surface number density from the optically selected sample can be explained by the radio 
selection as well as the reasons mentioned above. In any case, we expect that our sample is 
close to complete at the NVSS radio flux limit only for quasars at R < 16 and z < 0.4. 

Objects which are not classified as quasars are mostly Galactic objects, but some have 
featureless spectra similar to BL Lac objects. These objects will be discussed in a future 
work. 
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5. CONCLUSION 

By combining the radio and NIR quasar selection algorithm, we have derived an algo- 
rithm which is efficient for selecting quasars with a minimal contamination rate. We applied 
this algorithm to search for quasars in low Galactic latitude regions which are crowded with 
stars. From the first two observing runs in 2006, we discovered 40 bright quasars/AGNs out 
of 88 candidates. The newly discovered quasars/AGNs will be useful for various astronomical 
studies such as the distribution of Galactic matter and the check of reddening maps, and 
the quasar host galaxy study using a nearby bright star as a guide star. 
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gram of the Korea Science and Engineering Foundation. This work uses a CCD camera 
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Table I. Properties of the newly discovered quasars/AGNs. 
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Table 1 — Continued 
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SNUQSO Name 
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( 8 ) _ ( 10 )iJ, B, and i-ST-band apparent magnitudes, not corrected for Galactic extinction. 
( n ) J — K color, not corrected for Galactic extinction. 
' 12 ^i?-band absolute magnitude, corrected for Galactic extinction 
< 13 )Flux at 1.4 GHz from NVSS. 
< 14 ) Reddening parameter, E(B - V) 

(IB), (16) Extinction at B and R-bands in magnitude unit, respectively. B-band extinction is calculated using Ag = 3.91 E(B — V) for USNO magnitudes, 
and Ab = 4.3 E(B — V) for the Maidanak data. For Ar, we used Ar = 2.93 E(B — V) for both cases, since the difference in the extinction correction arising 
from a small discrepancy in conversion factor is less than 0.01 mag between USNO-II R magnitude and the Maidanak R magnitude. 

( lr ) Comments. 'X' stands for X-ray source, and 'GB' for Gamma-ray Blazar Candidate. 
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Fig. 1. — Distribution of quasars and AGNs in the sky, plotted in Galactic coordinates. 
Small dots represent the quasars/AGNs brighter than V < 23 mag discovered to date, while 
crosses represent quasars/AGNs brighter than V = 17 mag. Note a gap in the distribution 
between the northern and southern Galactic hemispheres. The newly discovered quasars are 
marked with squares (from 2006 June run) and triangles (from 2006 December run). [High 
quality image is available on electronic journal] 
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Fig. 2. — Optical spectra of the newly discovered quasars, plotted in the order of increasing 
E(B — V) values. The E(B — V) values increase from top to bottom, left to right. Broad 
Balmer emission lines and [O III] lines are visible for most of the quasars except for those at 
the location where the Galactic extinction is high. 
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Fig. 3. — .R-band apparent magnitude of the newly discovered quasars before (squares and 
solid line) and after (circles and dashed line) the extinction correction. From this figure and 
the analysis of the surface number density (see text), we find that our sample has an effective 
optical magnitude cut of R < 16 mag. 



